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ABSTRACT 

We investigate the consequences of a non-negligible baryon fraction for models of 
structure formation in Cold Dark Matter dominated cosmologies, emphasizing in par- 
ticular the existence of oscillations in the present-day matter power spectrum. These 
oscillations are the remnants of acoustic oscillations in the photon-baryon fluid before 
last scattering, for which evidence from measurements of the Cosmic Microwave Back- 
ground anisotropy is mounting. For acceptable values of the cosmological and baryon 
densities, the oscillations modulate the power by up to ~ 10%, with a 'period' in spa- 
tial wavenumber which is close to Ak ~ 0.05 Mpc" . We study the effects of nonlinear 
evolution on these features, and show that they are erased for k > 0.2 h Mpc -1 . At 
larger scales, the features evolve as expected from second-order perturbation theory: 
the visibility of the oscillations is affected only weakly by nonlinear evolution. No 
realistic CDM parameter combination is able to account for the claimed feature at 
k ~ O.f /i _1 Mpc in the APM power spectrum, or the excess power at 100 ft. -1 Mpc 
wavelengths quoted by several recent surveys. Thus baryonic oscillations are not pre- 
dicted to dominate existing measurements of clustering. We examine several effects 
which may mask the features which are predicted, and conclude that future galaxy 
surveys may be able to detect the oscillatory features in the power spectrum pro- 
vided baryons comprise > 15% of the total density, but that it will be a technically 
challenging achievement. 
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1 INTRODUCTION 

It is commonly assumed that the mass of the Universe 
is dominated by Cold Dark Matter (CDM) and that the 
baryons had a negligible influence on the development of 
large-scale structure. Recent estimates of the ratio of the 
baryon density to dark matter density suggest that this view 



reionization, in which case they could only be revealed by 
large galaxy surveys. 

The baryon density as determined by Big Bang Nu- 
cleosynthesis (BBN) has been o n the rise in recent years 
(see discussion in White et al. ( L996| )). At the beginning 
of this decade, the baryon density was estimated to have a 
95% confidence range of Q, n h 2 _ 0.0125 ± .0025 flWalkei 



may be incorrect. Because of the coupling between baryons 
and the Cosmic Microwave Background (CMB) radiation 
during the recombination epoch, additional structure in the 
power spectrum will develop compared with a standard 
CDM scenario in which the power spectrum contains only 
one feature, at the horizon scale of matter-radiation equal- 
ity. We discuss the detectability of such higher order effects 
in the context of currently favored models for structure for- 
mation. We show that they will leave a significant imprint 
on the matter power spectrum that will be measurable with 
forthcoming galaxy redshift surveys for a wide range of pa- 
rameters, including the possibility of a direct detection of 
acoustic oscillations in the galaxy power spectrum. Unlike 
the CMB, the imprint of the oscillations on the matter power 
spectrum will survive even if the universe underwent early 



ct al. 1991 Smith, Kawano fc Malaney 1993 ), where Qb is 
the baryon density in units of the critical density, and h is 
the present Hubble parameter in units of 100 km s _1 Mpc - . 
This range has recently broadened and shifted toward higher 
values. For example, the measurement of the pri mord ial deu- 
terium abundance by Tytler, Fan & Buries ( [t996| ) yields 
n B h 2 = 0.024 ± 0.002 ± 0.002 ± 0.001, with the la uncer- 
tainties being statistical, systematic, and theoretical respec- 
tively. While extragalactic assessments of Ob have crept up- 
wards, the measured abundanc e of deuterium in the ISM 
sets a ceiling of Qsh 2 < 0.031 (Linsky et al. 1995), within 



the assumptions of the standard BBN paradigm. It may be 
possible to relax the nucleosynthesis bound on Qb with new 
particle physics, an example being the decaying tau neutrino 
proposal of Gyuk & Turner (1994), or by allowing inhomo- 
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geneities in the baryon-to-photon ratio (see e.g., Mathews, 
Kajino & Orito 1996 and references therein). However, in the 
latter case, even with relaxed constraints on the primordial 
Li abundance one finds n B h 2 < 0.0325. 

Support for a hi gh baryon fraction i s also provided by 



The recent compilation by White & Fabian (199E 



1995;. 

199q) 

gives 
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clus ters of galaxies ( White et al. 1993h ; White 
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(95% confidence) with comparable l ower limits quoted by 
other authors: Steigman & Felten (1995) find Q,b/£Io > 
0.2 (/i/0.5)" 3/2 , while Evrard, Metzler & Navarro fll996| ) es- 
timate fi B /fto > 0.11 (/i/0.5)" 3/2 . Here fi = ^cdm + fi B 
is the total non-relativistic energy density relative to the 
critical density. Given the above discussion we consider 
0.01 < n B h 2 < 0.03 and 0.1 < n h 2 < 0.25 as a fair range 
of parameters and we shall work within this range except in 
Section |^ where we shall relax these assumptions. 

For some time it has been recognized that a universe 
with a substantial baryon fraction will result in signifi- 
cant features in the matter power spectrum. Historically, 
Sakharov predicted oscillations in the matter power spec- 
trum in a cold universe with no radiation component, build- 
ing on earlier work of Lifshitz. The oscillations were an im- 
print of sound waves as in the modern context, but the 
restoring pressure was due to degenerate electron pressure 
at high densities, not the CMB photons. The first calcula- 
tions to emphasize the oscillations in the context of hot big 
bang CDM models with isocurvature and adiabatic fluctu- 
ations concentra ted on mode l s with very low Op ~ 0.1 and 
high Q B ~ flo (fctekel 1984| |Bardccii, Bond fc Efstathiou 



1987]; |P|eebles 1987a|; |Peebles 1987b|; [Blumcnthal ct al. 1988f 



Sugiyama. 1995| ; [Hu fc Su" giyama 1996| ). The features are ex- 
tremely prominent in this case - but the parameter values 
involved in these models lie well outside of the currently 
preferred range. 

In this paper, we investigate more realistic values of 
the parameters where the features are correspondingly more 
difficult to detect. A study of the CDM + baryon transfer 
function over a wide range of the SIo-^b parameter space 
has recently been p erform ed by Eisenstein & Hu (199£) 



and Eisenstein et al. (1998), and the detec tability of oscilla- 



tions ha s been studied by Tegmark (1997) and Goldberg & 



Strauss (1998). Where our results overlap they are in good 
agreement; our work concentrates on the role of nonlinear 
growth in modifying the form of the oscillations, which has 
not previously been explored. 

The outline of the paper is as follows. In Section ^ 
we discuss the predictions for the oscillations in the mat- 
ter power spectrum from linear theory (some of the relevant 
technical details are in Appendi xjA| ) . The effects of non- 
linearities are included in Section Musing 2nd order pertur- 
bation theory and N-body simulations (our implementation 
of a PM code is discussed in Appendix ^) . In Section ^ we 
discuss the existing evidence for large scale features in the 
power spectrum, and conclude that baryonic features are 
unlikely to account for the existing power-spectrum data. In 
Section ^ we discuss the measurability of these features in 
forthcoming redshift surveys and in larger samples selected 
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Figure 1. A comparison of the power spectra for the CMB and 
LSS, for a model with C = 0.4, A = 0, h = 0.65 and C B = 0.045. 
The upper panel is contribution per ln(fc) to the total variance 
in the CMB (scaled to unity at the smallest wavenumbers) . The 
lower panel is the power per ln(fc) in the matter, i.e. A 2 (fc), with 
the trend taken out using the BBKS fitting function (see text). 
The vertical dotted lines mark the positions of the peaks in the 
CMB power spectrum, which are 90° out of phase with the cor- 
responding oscillations in the LSS power spectrum. 



via photometric redshifts. This section includes a discussion 
of expected error bars, plus systematic effects which need 
to be considered, some of the details of which are outlined 
in Appendices |c] and |d| We summarize our conclusions in 
Section |(| 



2 LINEAR THEORY 

There is considerable evidence to suggest that the matter 
density of the universe is less than critical (flo < 1). We have 
discussed above the growing evidence favouring a baryon 
fraction larger than has been assumed until recently. This 
has motivated us to look in detail at a variety of CDM mod- 
els, including those with flo < 1, and to examine cases with 
a baryon fraction consistent with current constraints. For 
some of this regime the gravitational potential is not totally 
dominated by the CDM and features in the matter power 
spectrum 



A L (fc) 



fc 3 P(fc) 
2tt 2 



5 »(s) 



3+n 



T 2 (k) 



(2) 



will be generated. Here 8u is the fractional density contrast 
(8p/p)k at horizon crossing fc = Hq/c, and T(k) is the trans- 
fer function. Specifically, the linear power spectrum A L (fc) 
contains a series of peaks of small amplitude which arise 
from acoustic oscillations in the baryon-photon fluid prior 
to recombination, the same source as the much larger peaks 
in the CMB power spectrum, for which there is growing ev- 
idence. 

We show in Fig. [j] the radiation and matter power 
spectra for a CDM model with Qq — 0.4, h = 0.65 and 
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QbIi = 0.02, calculated by numerical evolution of the cou- 
pled Einstein, Boltzmann and fluid equations. Note that the 
radiation power spectrum Ay is the contribution per In fc to 
the total variance in the CMB; it is closely related to the 
angular power spectrum £ 2 Ci, but the two should not be 
confused. The top panel shows the radiation power spec- 
trum, which exhibits a clear series of peaks. These peaks 
are modes which are density maxima and minima of the os- 
cillations of the photon-baryon fluid at recombination (see 
Appendix [X| for more details). The bottom panel shows the 
square of the matter transfer function T(k) for this same 
model, with the gross features (the bend from T(k) ~ 1 as 
k — *• to T(k) ~ k~ 2 at k — > oo) removed by dividing by th e 
fitting function for T(k) provided by Bardeen et al. (198C) 
(hereafter BBKS) 



T(k) = 



In (1 + 2.34?) 



2.34<7 



[1 



3.89<? + (16.1q) 2 + 
(5.46g) 3 + (6.71#r 1/4 , 



(3) 



with q = k/hT* and T* a parameter. The BBKS fitting func- 
tion provides a reasonable fit to T(k) over the range of scales 
plotted where T(k) changes by 2 to 3 orders of magnitude. 
It is also in wide use and provides a practical reference stan- 
dard. By taking the ratio of our results to the BBKS results 
we can focus on finer scale features in the spectrum, with the 
prominent trends due to matter-radiation equality removed. 

We set the parameter T* in the BBK S fitting function 
using the prescription of Sugiyama ( 1995j ) : 

r* = Q h cxp[-n B (l + V2h/Qo)]. (4) 

We have labeled this parameter F* rather than T to empha- 
size that it governs the shape of the transfer function alone, 
and not the shape of A 2 (fc) if the initial spectrum is not 
scale- invariant. A more i nvolved expression for T* is given 
by Hu & Sugiyama ( |l996| ) (Eqs. D29, E12), but we shall use 
the simple scaling; it has received widespread use elsewhere, 
and it is informative to examine its accuracy. The BBKS 
fitting function itself is accurate only to about 5%, even in 
the limit Qb — > 0, and this accounts for some of the dis- 
crepancy shown in Fig. |l| Finally we should emphasize that 
Fig. |l| does not show the full effect of including baryons in 
the calculation of the power spectrum. Compared to a zero 
baryon model with the same Qo, a model with baryons has 
less power on all scales smaller than k ~ 10 -2 /iMpc -1 . This 
suppression of power on small scales is partly accounted for 
by the rescaling of T in Eq. ^. We reiterate that we have 
divided by the BBKS form purely in order to remove the 
trend and focus on the finer features in the power spectrum. 

Fig. ^ shows that on scales k ~ 0.1 /iMpc" 1 , there is a 
suppression of power relative to the BBKS fitting function. 
Partly this is due to inaccuracies in the BBKS fit, but there 
is also a physical effect at work: the baryons are supported by 
photon pressure at early times and their fluctuations cannot 
grow in amplitude until they are released from the photons. 
On small scales the peaks in the matter and radiati on power 
spectra are out of phase due to velocity ov e rshoot ( gunyacv 



& Zel'd ovich 197C ; Press fc Vishniac 1980; Hu fc Sugiyama 



199f). The growing (decaying) mode of the perturbations 



projects primarily onto the velocity (density) of the pertur- 
bation at high k. At larger scales the phases of the peaks 
become comparable, as the growing mode is sourced more 
by the density. Further discussion is given in Appendix |X| 
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Figure 2. The position (in k space in Mpc -1 ) of the second 
maximum in T(k) compared to BBKS, as a function of Clfjh 2 and 
!)b/i 2 . Here Oo = ^CDM + £^B- For measurements of distances in 
h Mpc -1 the values of the contours need to be divided by h. Apart 
from this scaling, the position of the peak is almost independent 
of cosmological parameters. 



A thorough discussion of the physics in linear theory and an 
analyti c fitti ng function for T(k) are presented by Eisenstein 
& Hu (1998) an d the application to 100- Mpc clustering by 
Eisenstein et al. (1995). Where our results overlap, they are 
in good agreement. 

The positions of the peaks in P(k) are very insensitive 
to the cosmological parameters over the range of interest. In 
adiabatic CDM models the second peak is the most likely 
to be measured first, since it occurs at a more accessible 
scale than the first peak. We show the position (in Mpc -1 ) 
of the second peak in adiabatic CDM models, as a func- 
tion of floh 2 and Q,Bh 2 in Fig. ^| The transfer functions 
for a grid of 35 models in total were calculated by numeri- 
cal integration of the coupled Einstein, fluid and Boltzmann 
equations from before equality through to the present. In 
Fig. |^ we show T{k) for a set of these models. On the ab- 
solute scale of the power spectrum, the features are small. 
Relative to the BBKS transfer function, however, the fea- 
tures are prominent. The peaks in Fig. ^Jwere found as the 
maxima and minima of T(k)/TBBKs{k) with an additional 
trend removed to account for the dip near fc ~ 0.1 /iMpc - , 
as is visible in Fig. ^. Notice that the position varies by only 
10% for reasonable variations in these parameters, and can 
thus be very reliably predicted by models. For isocurvature 
models the peaks are shifted by a factor of approximately 
1.5 to higher k. The other peaks form an almost harmonic 
series, as shown in Fig. ^ and discussed in Appendix [X]. 



3 BEYOND LINEAR THEORY 

The power spectra shown in Figs. |l| and ^ are computed 
using linear perturbation theory, as discussed in Appendix 
|X|. However for COBE normalized models, even the first 
few features are in the trans-linear regime (A^ ~ 0.1-1), 



© 1998 RAS, MNRAS 000, 000-000 



4 A.Meiksin et al. 




0.01 0.1 1 0.01 0.1 1 



k (h/Mpc) k (h/Mpc) 

Figure 3. The transfer function T(k) for Ho = 0.3, 0.4 and h = 0.6, 0.65 and 0.7. Note that the features in T(fc) are quite small on the 
scale of the variation of T(k), but prominent when the trend is removed, as shown in Fig. kl 
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Figure 4. The linear theory matter power spectrum with the trend removed by dividing by the best fitting BBKS analytic fit to T(k). 
Higher gives larger oscillations. 



so it is far from clear that linear theory will apply to these 
oscillations. To compute the full nonlinear power spectrum 
A 2 , we have gone to second order in perturbation theory 
and performed N-body simulations as described below. Our 
conclusion is that the general trend of the non-linear effects 
is to suppress the peaks in the non-linear regime (Figs. ^, Q 
§and§). 

3.1 Normalization 

Since the persistence of the peaks depends on whether they 
lie above or below the scale of non-linearity, the absolute nor- 
malization of the models is important. The COBE normal- 



ization distinguishes between open and vacuum-dominated 
models. In the absence of tensor contributions, the horizon 
crossing amplitude is 

Su = 1-95 X 10 -5n-°-35-0.191n no -0.17(n-l) x 

exp[-(n- 1) - 0.14(rt- l) 2 ] 
for open models and 



(•») 



<5W = 1-94 x 10" 5 a 



-0.785-0.05 lnf! 



exp[-0.95(?i - 1) - 0.169(ra - l) 2 ] 



(6) 



for flat models flBunn fc White 1997| ). Thus if the models 
are COBE normalized we would expect the open models to 
show more pronounced oscillations. However one can also 
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Figure 5. The values of primordial spectral index required in 
order for CDM models to satisfy the constraints of correct nor- 
malization on both COBE and cluster scales, for the indicated 
values of h. This figure assumes no tensor contribution to the 
CMB anisotropy. 



Model 


Oo 


n B h 2 


h 


n 


r* 


CT8 


tCDM 


1.0 


0.0375 


0.50 


0.748 


0.370 


0.55 


OCDM 


0.4 


0.0300 


0.65 


1.258 


0.198 


0.92 


ACDM 


0.4 


0.0300 


0.65 


1.030 


0.198 


0.92 


EHSS 


1.0 


0.1440 


0.60 


0.889 


0.259 


0.55 



Table 1. Parameters for the cosmological models considered in 
the text. The column labeled T* is the parameter T of the BBKS 
transfer function which best fits the numer cal t ransfer function 
of the model, as determined by Sugiyama (1995). These are not 
the values of T which would be inferred from A 2 (k) since n =i 1. 
The model EHSS is the high-Qe model suggested by Eisenstein 
et al. (1998), with n reduced in order to fit the cluster abundance. 



normalize the models on small scales (k ~ 0.2 ftMpc -1 ) us- 
ing the abundance of rich clusters (see e.g. Pen 1998 and 
references therein) . In this c ase the open and A mode l s have 
very si milar normalizations ( Eke, Cole fc Frenk 1996 ; Viana 
& Liddl e 1996| ; [Eke, Cole, Frenk fc Henry 1998j ; [Viana fc Lid- 
dle 199 1). We adopt the following simple fitting formulae 
which is near to the mean opinion on this issue: 

as =0.55^ 56 (7) 



(White, Efstathiou fc Frenk 1993a). In what follows, we will 



generally force the models to fit both COBE and the abun- 
dance of rich clusters by adjusting the primordial spectral 
index n appropriately. Fig. H shows the values of n required 
in order for CDM models to satisfy these constraints for var- 
ious values of flo, and h. This was calculated usi ng th e 
approximate transfer function of Eisenstein & Hu (1998), 
and neglecting any tensor contribution to CMB anisotropies. 
It is apparent that, for reasonable baryon content, £7o — 1 
models generally require a substantia l tilt towards n < 1 , 
whereas low-fio models need n > 1 (White & Silk 1996). 
The smallest degree of tilt is required in the case of ACDM 
and a low baryon fraction. 
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Figure 6. The effect of non-linear evolution on the oscillations 
in the power spectrum. We show (solid line) the ratio of the 
line ar pow er spectrum for the model suggested by Eisenstein et 
al. (1998) with O = 1, h = 0.6, Cl B = 0.4 and n = 0.889 (we 
have lowered n to provide a better fit to the cluster abundance) 
to a BBKS model. The BBKS transfer function has T* = 0.259 
and the model is COBE normalized. Note that the BBKS transfer 
function with T* as provided by Sugiyama (1995) does not provide 
a good fit to the numerical T(k). At small scales a much better fit 
is provided by T* = 0.265. The dashed line shows the ratio includ- 
ing 2nd order corrections, the small arrow near k ~ 0.4/iMpc -1 
indicates where A 2 (fc) = 1. The open squares are the results from 
N-body simulations with the CDM spectrum in real space, solid 
squares the results in redshift space. Note that the oscillations for 
k > 0.1?tMpc -1 are washed out. 



3.2 Perturbation Theory 

Since the features we are interested in occur in the regime 
where A 2 (fc) ~ 1 we first computed the correction to the 
l inear theory resul t using second o rder perturbation theory 



(^brszkiewicz_l^M 



., [Vishniac 1983| ; |Makino et al. 1992| ; [Iain 
fctBert^chmgeTT 994| ) . The basic method involves evaluating 



two convolution integrals numerically, 



A? a) (fc)=A£(ft) 



+ / dlnqAUq) 
hi 



A£(|k-q|)Jir 1 (k,q) 



+Al(k) / d\nqAi(q)K 2 (k,q), 
Jo 



(8) 



where A 2 2 j(fe) is the power spectrum including second or- 
der corrections, and K\(k., q) and K2(k,q) are integration 
kernels. These are given by 



*x(k,q) = 



k 



196 <? 4 |k-q| 7 

,2,2 , ~,3. 



(3g 2 fc 2 + 7qk :, fi - Wq 2 k 2 ^ 2 ) 2 
where /i = k • q/(fcg), and 

252 q 2 



(9) 



K 2 (k,q) 



k 2 a 2 o 4 

V^ 79 + 50 fr~ 21 F 
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Figure 7. As in Fig. |6J but for a model with Qo = 1, h = 0.50, 
C B = 0.15 and n = 0.748. The BBKS transfer function has T* = 
0.37 and the model is COBE normalized. 



Figure 8. As in Fig. g but for a model with Qo = 0.4, Qjy = 
0.6, h = 0.65, n B h 2 = 0.03 and n = 1.030. The BBKS transfer 
function has T* = 0.198 and the model is COBE normalized. 
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The perturbation expansion formally breaks down for large 
A 2 (fc). In principle the integration over q should be termi- 
nated at some maximum wavenumber; in practice, the re- 
sults are not very sensitive to the upper limit. We note that, 
for an Einstein-deSitter universe, the second order correc- 
tion to the power spectrum evolves as a 4 , while the power 
spectrum in linear order evolves as a 2 , where a is the ex- 
pansion factor. We adopt the same expressions for the open 
models below, since the additional corrections to the second 



order terms appear small ( Catelan et al. 1995 ). 

One may also compare the second-order calculation 
with the fully n onlin ear scaling formula developed by Pea- 
cock & Dodds (1996). Extension of this formalism to treat 
oscillatory spectra is difficult since there is no unique defi- 
nition of a local spectral index. A naive application of the 
results, using the spectral index from the best-fitting BBKS 
spectrum and the full oscillatory power spectrum in A 2 , re- 
sults in a very poor match to the exact non-linear spectrum. 
We have not pursued this matter further. 



3.3 N-body models 

In order to find the exact behaviour of the spectrum, we 
ran a series of N-body simulations using two P 3 M codes. 
The first was an adaptation of Hugh Couchman's AP 3 M 
code (Couchman 1991; see also Peacock & Dodds 1996) and 
the second was the Hy dra N-body, hydrodynamics code, ru n 
in collisionless mode (Couchman, Thomas & Pearce 1995). 
The agreement in A 2 (fc) between both of these codes and 
an independent PM code (see Appendix |§|) was very good, 
indicating that the numerical results were stable. Note that 
these codes treat all matter as cold, whether or not some 
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Figure 9. As in Fig. g but for a model with Qo = 0.4, Qjy = 
0, h = 0.65, fl B h 2 = 0.03 and n = 1.258. The BBKS transfer 
function has T* = 0.198 and the model is COBE normalized. 



fraction is baryonic or HDM. Such differences have a strong 
effect on small scales, but for k < 1/iMpc -1 are negligible. 

Having determined that the PM code adequately calcu- 
lated the non-linear power spectrum on the scales of interest, 
we used its superior speed to allow the averaging over many 
realizations. For each model the PM code was run many 
times (typically 10 3 in 64 3 mode for the low-fc modes and 
fewxlO 2 in 128 3 mode for the high-fc modes) with a start- 
ing redshift of 1 + 2 = 20, and different Gaussian initial 
conditions. The resulting determinations of A 2 (fc) were av- 
eraged over the runs. An alternative procedure designed to 
reduce scatter at low-fc, allowing random phases but with 
fixed amplitudes (Peacock & Dodds 1996), was found not 
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to be sufficiently accurate for our purposes. With a careful 
treatment of binning effects, we believe our results for the 
power spectrum are accurate to 1%. 



3.4 Real— space power spectra 

The results of these procedures are shown in Figs. ^J-^)| The 
high-fc oscillations in A 2 (fc) are suppressed once second or- 
der effects are included. There is an additional reduction of 
power on somewhat longer scales (k < 0.1 ft Mpc" 1 ). Sec- 
ond order perturbation theory provides a reasonable, but 
not highly accurate, description of the non-linear effects. It 
successfully predicts the scale of complete suppression of the 
peaks, but breaks down at A 2 (fc) < 1. 

The clear message of these figures is that only one un- 
ambiguous baryonic feature in the spectrum is expected. The 
first baryonic peak relative to BBKS is hardly altered by 
nonlinear evolution, but it is not clear that it is detectable. 
As we discuss below, the likely future errors on the power 
spectrum at these very large scales will not allow such a 
10% feature to be picked out. In any case, it is very broad in 
fc-space extent, and will depend rather sensitively on which 
BBKS model is used as a reference. In contrast, the sec- 
ond peak is relatively narrow, being confined to roughly a 
factor 1.3 in k. This feature is usually near the maximum 
wavenumber at which it can survive nonlinear damping, and 
the observational errors are much smaller. Finally, we have 
shown in Section 2 that the position of this peak is quite ro- 
bust once h is determined. Future efforts to prove the impor- 
tance of baryons in determining the form of the fluctuation 
should therefore concentrate on this signature. 

The results so far have been entirely in real space. How- 
ever, the only statistics which are independent of redshift- 
space effects are projected statistics, e.g. the projected cor- 
relation function: 



£[(r 2 +x 2 ) 1/2 ] dx, 



(11) 



or the angular correlation function w(0). Integration over a 
wide selection function damps the visibility of oscillations. 
With typical survey selection functions the features are sup- 
pressed to an unmeasurable level. To avoid this one must 
use sub-samples with a narrow distribution in depth. We 
will discuss the Fourier space analog of H(r) calculated on 
narrow shells selected using photometric redshifts in Section 



3.5 Redshift— Space Power— Spectra 

The redshift-space power spectrum is distorted from the 
real-space power spectrum by peculiar ve locities. On large 
scales this leads to an increase in power (Kaiser 1987) be- 
cause objects stream towards overdensities. On short scales 
virialized motions suppress clustering in redshift space. 

We have computed the power spectrum in redshift 
space, A 2 (fc), from the N-body codes. For simplicity, an in- 
finitely distant observer is assumed; i.e., the redshift dis- 
tortions are along a single Cartesian axis of the box. The 
results are in reasonable agreement with analytical mod- 
els of the re dshif t space distortion which involve convolving 
the Kaiser (1987) result with an exponential on small scales 



( |Cole, Fisher fc Weinberg 199*51 ). However, the transition 
from enhancement to suppression of A 2 (fc) starts to set in 
nearfc~ 0.1-0.2 h Mpc" 1 , and by k ~ 0.3-0.5 h Mpc" 1 the 
suppression is severe. Since the position of the features we 
are proposing to observe are near k ~ 0.1 ftMpc -1 , an inter- 
pretation of the power spectrum in redshift space requires a 
careful modeling of the redshift-space distortions. 

We show A 2 (fc) as the solid squares in Figs. [], (?], [| and 
Note that the appearance of the baryonic features are al- 
tered by the redshift space distortions, but persist in a rec- 
ognizable form. As in the real-space results, the only feature 
that will be robust with respect to how the smooth under- 
lying spectrum is defined will be the second peak, which is 
expected to give a roughly 10% boost to the power in a range 
of k of about a factor 1.3 centred on k ~ 0.055 Mpc -1 . The 
detection of this signature will be an observational challenge 
but, as we discuss below, it may be feasible. 



4 FITTING CLUSTERING DATA 

So far our discussion has been confined to a firm prediction 
for the CDM family of models: for large Qb and low Slo, 
we expect to see oscillations in the matter power spectrum 
arising from the same mechanism (acoustic waves before de- 
coupling) which produced the peaks in the CMB spectrum, 
for which evidence is mounting. We next consider the pos- 
sibility that tentative evidence for the oscillations has been 
found in recent surveys. 

There are several surveys which have quoted evidence 
for excess po wer on 100- Mpc scales. The first was the Broad- 



hurst et al. (1990) pencil beam surveys, which quoted an 



excess of power on 128 h 1 Mpc wavelengths . The 2D anal- 



ysis of the LCRS survey (Landy et al. 1996) also reported 



a statistically significant bump in the power spectrum at 
k ~ 0.06 /iMpc" 1 (b ut see Bromle y & Press 1998). Re- 
cently Einasto et al. (1997a; 1997b) have claimed to find 



120 h~ x Mpc periodicity in the 3D distribution of superclus- 
ters. Evidence for a distinct kind of excess power on large 
scales is the broad bump detected around k ~ 0.1 feMpc -1 



in the APM and IRAS real- 


space data ([Saunders, Rowan- 


Robinson & Lawrence 1992; 


Maddox, Efstathiou & Suther- 


land 1996; Peacock 1997). 



The intriguing possibility is that these surveys may 
be indicating the existence of acoustic oscillations in the 
underlying matter power spectrum. From our earlier dis- 
cussion it appears that one cannot explain these features 
from the linear-theory mass power spectrum of an adiabatic 
CDM model with Ob constrained by standard BBN. What 
would be required to explain such features? This question 
has be en ad dressed recently in linear theory by Eisenstein 
et al. (1998). They identify two models which might fit the 
data, one with low fio and one with high fio- As expected 
from the discussion in Section [|, the parameters required for 
either of these models do not appear credible from other con- 
siderations (e.g., BBN constraints on the baryon fraction). 
However these models provide, by construction, features in 
the places for which there is (controversial) evidence as dis- 
cussed above. 

To see how the features are affected by non-linearities, 
we have run N-body simulations of the high-fi) model 
(fio = 1, fiB = 0.4, h = 0.6, n = 0.889). The result is shown 
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in Fig. |(| and differs somewhat from the other models we 
have investigated. Because of the extreme parameter choices 
in this case, the baryonic features are at somewhat larger fc. 
This aspect allows the model to place a broad feature of am- 
plitude roughly a factor 1.8 at the desired fc ~ 0.05 /iMpc -1 
- but the second and subsequent peaks are then at suf- 
ficiently high k that they are destroyed by nonlinearities. 
There has to be serious doubt over whether such a rela- 
tively broad 3D feature co uld cause a significant effect in 
lower-dimensional surveys (Kaiser & Peacock 1991), but in 
any case we have shown that smaller-scale harmonics will 
not survive in order to give a test of such a hypothesis. 

An effect that may alter the relative heights of the peaks 
is the relative bias between different galaxy populations that 
have different redshift distributions. By treating the two 
populations as a single population, a scale dependent bias 
in the power spectrum will be introduced in the measured 
2D power spectrum, because the angular scale correspond- 
ing to a given 3D clustering strength will differ for the two 
populations. For example, a luminosity-dependent cluster- 
ing amplitude could lead to such an effect. 

The standard prescription for calculating the power 
spectrum from a survey implicitly assumes that the clus- 
t ering is independent of the lumino s ity. There is evidence 



(Davis et al. 1988} Lin et al. 1996; Bromley et al. 1998 



Willmer et al. 199q ) that this is not the case. We have tested 
the effects of two populations with differing bias factors and 
selection functions on the calculation of the power spectrum, 
presuming the two populations are treated as a single pop- 
ulation. We find that in all cases the effect is weak. The 
only case in which a large effect can be seen is if one of the 
populations has a tight selection function in z, though we 
regard this case as artificial. For realistic selection functions 
the scale dependent bias introduced by two populations with 
61/62 "S 2 is at the 1% level. We give details in Appendix |c| 
While the distribution of power in a single mode is ex- 
ponential even in the presence of strong nonlinearities ( |Fan| 
& Barceen 1995), in a finite- volume survey nonlinear con- 
tributions may increase the dispersion i n the power est i- 
mate above the exponential expectation (Amendola 1994). 
This may be viewed as a consequence of the finite num- 
ber of statistically indep endent spatial cells in any given 
survey. Amendola ( |1994| ) showed that this may be a par- 
ticularly strong effect in pencil-beam surveys. An extension 
of the analysis to higher dimensions shows that the non- 
exponential terms decrease in size inversely with the num- 
ber of independent cells within the survey. The effect is also 
in part a consequence of the small number of modes avail- 
able to measure power at low k. Since modes are angle- 
averaged within frequency bands in multidimensional sur- 
veys, a large number of independent modes will contribute 
to band-averaged power spectrum estimates even on large 
scales. This ensures that the power estimates will follow a 
X 2 distribution with a large number of degrees of freedom. In 
the limit that the number of independent modes m — > 00, 
[P(fc)/(P(fc))] 1/2 , where P(fc) is the band-averaged power 
near k and (P(fc)} is the expected power, tends to a nor- 
mal Gaussian deviate wit h unit mean and variance l/4m 
(Kenney & Keeping 1959). Thus it would appear that any 
features in the power spectrum on the scales of interest are 
not likely to be large deviations from the mean power on 
those scales. 



5 FUTURE OBSERVATIONAL PROSPECTS 

Since the features that we are proposing to measure are 
only small perturbations to the power spectrum, they are 
difficult to see in current surveys (see Section ^). However 
both the AAT-2dF and the Sloan Digital Sky Survey (SDSS) 
will soon enhance our knowledge of the 3D geometry of the 
local universe. We focus here on the SDSS. 



5.1 3D Surveys 

The SDSS Northern Polar Cap redshift survey will contain 
10 6 galaxies with a median distance of ~ 350 /i -1 Mpc. An 
additional ~ 5 x 10 7 galaxies will be measured in 5 bands 
in a photometric survey ~ 5 magnitudes deeper than the 
redshift catalog. It will be possible to measure the 2D power 
spectrum (or angular correlation function) using this deeper 
survey. It will additionally be possible to extract measure- 
ments of the 3D power spectrum from the larger catalog by 
the use of photometric redshifts. 

A critical issue for the detection of baryonic features is 
the expected precision of these power estimates. Feldman et 
al. (1994; FKP) studied this question in some detail, and 
we now summarize their main results (see also Tegmark et 
al. 1998) . Throughout it is assumed that on the scales of in- 
terest the modes remain approximately independent. While 
this is roughly true, in order to attain the high levels of 
precision of which the SDSS is capable this a ssumption may 
need to be improved (|Meiksin fc White 1998]). The FKP for- 



mula for the fractional variance in the power measured for a 
Gaussian field by averaging modes in some fc-space volume 
Vfe is their equation (2.3.2): 



(2tt) 3 / d 3 r n 4 w 4 [l + l/nP] : 



V k 



[J d?r 



(12) 



where n(r) is the mean density defined by the survey selec- 
tion, and w(r) is a weight function, which should be set to 
w = [1+nP] -1 for minimum variance. This formula assumes 
that the fc-space region is large in all directions compared to 
the window function resulting from transforming the survey 
geometry. This assumption brea ks dow n at the very largest 
scales; see Goldberg & Strauss (1998) for a more sophisti- 
cated analysis. Also, in the case where Vk is a shell in fc space, 
rather than an isolated region, the above formula should be 
multiplied by a factor 2 to allow for the fact that the density 
field is real (FKP equation 2.4.6). 

In the simple case of a large survey where shot noise is 
negligible, nP 3> 1, this expression becomes simply 



(^) 2 = [2]x(2n) s (V k Vr 1 , 



(13) 



where V is the survey volume. This reflects the fact that 
the density of states in k space is (2ty) 3 /V. The weight func- 
tion automatically cuts the survey off at the point where 
shot noise is starting to become important. We may choose, 
however, to analyze a volume-limited subsample for which 
the shot noise is negligible. 

FKP also give an expression (their 2.2.6) for the co- 
variance of the power estimates at different wavenumbers. If 
shot noise is negligible, this says that the two-point corre- 
lation function for fluctuations in power is just the Fourier 
transform of n 2 w 2 . However, for a shell whose width Afc is 
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large enough for the above expression for the variance in 
the shell-averaged power to apply, this correlation function 
is guaranteed to be small. The power estimates in adjacent 
bins will thus be effectively uncorrelated. 

As an example, the FKP effective volume for the 
Sloan survey is approximately 0.17 (/i _1 Gpc) 3 , taking P — 
5000 (/i~ 1 Mpc) 3 as a typical reference power, and u sing the 



selection function given in Gunn & Weinberg (1995). This is 
equivalent to uniformly sampling the volume out to 2 ~ 0.21; 
for an illustrative discussion, we henceforth assume that the 
Sloan survey will sample this volume in a way that is ef- 
fectively independent of shot noise (we do include the shot- 
noise component in the error bars, but it is a very small 
correction). For Q — 1 and the Sloan area of tv sr, the co- 
moving volume is equivalent to that of a sphere of radius 
approximately R = 340 h~ x Mpc. If the survey window was 
in fact perfectly isotropic, the power correlation function 
would be 



{SP(k)8P(k + Ak)) oc |/(-RAfc)| 



(14) 



where f(y) = 3(siny — y cos y)/y 3 . This has its first zero at 
Afc ~ 4.5/ R, suggesting that samples of the power at a spac- 
ing of > 0.013 h Mpc -1 will be uncorrelated. The transform 
of the real conical window function is in practice about 15 
per cent broader than would be calculated just on the ba- 
sis of the effective volume of the survey. Overall, this implies 
that we can treat estimates of power separated by more than 
about 0.015 ft Mpc -1 as being effectively completely uncor- 
related. Note that this is a rather more conservative crite- 
rion than t he A fc = 0.0043 h Mpc -1 assumed by Goldberg 
& Strauss (1998). According to our calculations, the power 



correlation coefficient at this separation should be around 
0.5. This may explain why they found slightly too small un- 
certainties in model fits to the power spectra when using 
FKP error bars. 

We show in Fig. [l(] the error bars which the SDSS 
Northern Polar Cap (NPC) redshift survey should be able to 
produce on the power spectrum. The error bars are for bins 
with a constant spacing of Afc = 0.015 h Mpc -1 in wavenum- 
ber. The above formula suggests a fractional power error of 



op 
P 



4tt 2 
Vfc 2 Afc 



= (fc/0.0039/iMpc~ 



(15) 



The NPC survey may be extended in several ways which 
will enhance its ability to measure A 2 (fc) at low fc. First will 
be by autocorrelation including the Southern Redshift Sur- 
vey, which provides several long baselines between the north 
and south to constrain the large-scale power. Secondly a 
deeper survey will be done using Bright Red Galaxies, which 
could reach to an effective depth of close to 1 Gpc. This will 
at least double the number of independent modes compared 
with the NPC survey, with a corresponding reduction of 
the errors by a factor of 2. In Fig. hoi we show representa- 
tive error bars corresponding to Afc = 0.0075 h Mpc -1 and 
ap/P = (fc/0. 00195 /iMpc -1 ) -1 . 

Because of the great depth to which QSO sources may 
be detected, it might be hoped that they could be used to 
obtain an even more precise measure of the power spec- 
trum. Unfortunately their sparse numbers would restrict 
any such measurement to being shot noise dominated in the 
forseable future. For example, Warren, Hewett, & Osmer 
(1994) find the proper density of QSOs in a sample slightly 



deeper (m or < 20) than either the SDSS or 2dF samples 
peaks at a value of uqso ~ 10 -3 ' 5 /i 3 Mpc -3 at z = 3.3 
(go = 0.5). A fiducial value for the power spectrum today 
of P = 5000/i~ 3 Mpc 3 corresponds to P w 300/i~ 3 Mpc 3 at 
z — 3.3, an order of magnitude smaller than the shot noise 
contribution n,Q SO « 3300. It is quite likely that the QSOs 
are a biased tracer of the underlying mass fluctuations. Still, 
the number of QSOs must be at least a factor 100/6qsq 
larger, where bQso is the bias factor, to detect the baryonic 
signature in the power spectrum. At this time it is unclear 
whether QSOs even exist in such numbers for plausible bias 
factors. 

Finally one could imagine using photometric redshifts 
for the full 5 x 10 7 galaxies in the photometric catalogue 
to compute A 2 (fc) at low fc where the effects of the redshift 
errors are less important. The larger number of galaxies al- 
lows one to probe deeper, and reduces the errors at low fc 
as in the BRG sample. How well this may be done in prac- 
tice will depend on the redshift distribution of the galax- 
ies. The distance errors from photom etric redshifts at low z 
are on the order of a ~ 50 hT 1 Mpc (Brunner et al. 1998). 



The power spectrum is strongly affected when fccr > 1, or 
fc > 0.02 h Mpc" 1 . On the other hand, if there is a suffi- 
cient number of galaxies to high redshifts, say z fa 1, then 
the decreased comoving differential path length with red- 
shift will increase the usable region in (comoving) fc-space 
up to fc « 0.04 ft Mpc -1 if errors as s mall as Sz — 0.02 ma y 
be achieved to these high redshifts (Connolly et al. 1995). 
This may be adequate for a detection of the primary, and 
largest, peak. Systematic errors are still potentially large in 
this range, but they may be sufficiently controllable. How to 
deal with the effect of evolution in these very deep surveys, 
however, needs further study. 

These questions aside, the error bars shown in Fig. |l^ 
represent a conservative estimate which may be able to be 
improved upon with more work. Clearly the ability of the 
SDSS to measure the oscillations in A 2 (fc) will depend on 
the underlying cosmological parameters. 



5.2 Degenerate models 

We quantify the ability of the SDSS to detect the presence 
of the acoustic oscillations by computing the expected \ 2 
between a power spectrum with oscillations and a corre- 
sponding limiting case with no baryons. In so doing, we as- 
sume that non-linear mode coupling can be neglected and 
the power estimates measured at the various sampling points 
are statistically independent of each other. Statistical inde- 
pendence is a stronger assumption than the result proved 
above, which is that sufficiently well-separated power esti- 
mates are uncorrelated. In practice, the difference appears 
to be small, and the reduced \ 2 values for the fit to the ex- 
act spectrum are very close to unity. We therefore neglect 
this technical distinction. Our neglect of non-linear mode 
coupling is slightly optimistic, but will not alter our results 
too much pr ovided we work at fc < 0.1 h Mpc -1 (Meiksin & 
White 1998| ). 

In Fig. |l(| we show a comparison between the fiducial 
ACDM model and a variant with n — 0.96, h — 0.56 and 
Qb ~ 0. The non-linear spectrum was again calculated from 
the N-body code, even thoug h the re are no oscillations to 
render the Peacock & Dodds (199£) procedure inapplicable. 
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timate on these scales, obtaining the low k measurements 
may be required. Alternatively, under the assumption that 
the dominant errors in both experiments are statistical, it 
is possible that a combined fit of the CMB anisotropy data 
from MAP and the SDSS data will break any degeneracy in 
the quasi-linear regime. This has been shown using a linear 
analysis by Eisenstein, Hu & Tegmark (1998). Note, how- 
ever, that the second harmonic in Figure h(j is invisible even 
to the BRG sample: over 0.045 < k < 0.1 h Mpc" 1 the BRG 
sample rejects the non-baryonic model with a confidence of 
only 90%. If the depth of the BRG is able to reach fully to 1 
Gpc, then the situation improves. With the increased num- 
ber of modes and precision, the BRG sample would then 
be able to distinguish the baryonic model from the non- 
baryonic model over the range 0.045 < k < 0.1 ft Mpc -1 
with a confidence of 99.93%. This optimistically presumes 
that evolution of BRGs over such a broad range in ages 
is well understood. It also neglects intrinsic correlations in 
the power spectr um, which are moderate ly large even near 
fc = O.lhMpc -1 flMeiksin fc White 1998|). 



Figure 10. A comparison of the fiducial ACDM model (solid 
line) with a variant of the model having no baryons (dashed 
line). The dotted error bars are estimates of la errors in uncor- 
rected bins which should be measurable with the SDSS North- 



Two effects need to be considered further before a case 
can be made for the detectability of baryonic features. The 
first is the question of whether bias alters their visibility. We 
have assumed linear bias, i.e. 5N g /N g = bSpm/pm, where N g 
denotes the galaxy counts, p m is the total matter density, 
and b is the bias factor. However, while semi-analyti c mod- 



eling has made some progress in u nderstanding bias ( Kauff- 



ern Poh r Cap survey assuming the underlying model is ACDM. 



mann, Nusser & Steinmetz 1997), it is not a priori clear 



The smaller solid error bars are for the Bright Red Galaxy sam- 
ple. The BRG sample will be able to detect the baryon sig- 
nature in the quasi-linear to linear regime only if the errors 
shown for k < 0.045 /iMpc - 1 are realized. The fiducial model 
has fi = 0.4 = 1 - fi A , tt B h 2 = 0.03, h = 0.65, and n = 1.030. 
The non-baryonic model has n = 0.96, h = 0.56 and Qb = 0. 
Both models are COBE normalized, and divided by the same 
BBKS model fit used in Fig. ^. (The small "oscillations" in the 
non-baryonic model at 0.03 < k < 0.06 h Mpc -1 are due to the 
limited accuracy of the non-linear calculation at low k, and are 
not physical.) 



We found that the non-linear scaling prescription was not al- 
ways accurate enough for our purposes, which may be an im- 
portant consideration when analyzing future high precision 
surveys. The model, which we fit by eye, has the large-scale 
clustering amplitude held fixed at the COBE central value 
(|Bunn fc White 1997]). We allowed the tilt and "shape" to 



vary to obtain the closest fit to our ACDM model. The data 
points and errors are those expected for the SDSS Bright 
Red Galaxy sample and the North Polar Cap sample, assum- 
ing ACDM is the true underlying model. The NPC sample 
is able to distinguish the two models at the 3cr level only 
for k > 0.33 /iMpc -1 . The greater depth of the BRG sam- 



whether the transformation between mass and light affects 
the relative visibility of the baryonic features. For the ACDM 
model we have carried out experiments with simple nonlin- 
ear modifications of the final density field, such as p — *• p B . 
This confirms that, on sufficiently large scales, linear bias is 
indeed a good approximation, even when the spectrum con- 
tains features (see also Scherrer & Wein berg 1998). Some 
results for p B "biasing" are shown in Fig. 0. Note that on 
large scales the bias is linear, while it becomes more scale 
dependent on small scales. Even for large enough B that 
the scale dependence would be strong, the second peak in 
the power spectrum is still present. These experiments re- 
main preliminary. A definitive answer to this question would 
require an understanding of the nature of galaxy bias. As 
a first step beyond our simple model one could compute 
the halo-halo correlation function, or include hydrodynam- 
ics to identify sites of galaxy formation. Neither of these 
approaches is feasible with our PM code. 

Secondly, when comparing power spectra measured at 
different redshifts, care must be taken to account for evo- 
lution. Within the biased galaxy formation paradigm the 
with z can be quite complicated and 



As has been 



pie per nits a clear distinction at the 3a level to be made 



b etweei|L the two models from the k < 0.045 faMpc ^mea 



surements. These measurements, however, may be difficult 
to achieve because of correlations in the photometry on these 
scales either due to intrinsic measurement error or as a re- 
sult of residual correlations in Galactic extinction. Confin- 
ing the comparison to larger k values requires points with 



evolution of Ag al 

non-monotonic, unlike the evolution of A„ 
pointed out r ecently by several authors (Brainerd fc Vil- 
_hmiscnT994 ; lOgawa, Roukema fe Yamashita 1997 ; Bagla 
1995 ; Prirnacket^iL 1998j ) the contribution to the clustering 



k > 0.1 j/tMpc" 1 , for which A J (fc) > 0.5, for the BRG sam- 



strength comes both from the underlying mass power spec- 
trum (which grows with z) and the bias of the object which 
is a function of its rarity at a given z (and typically decreases 
with z). Preliminary o bservational evidence for evolution has 
recently been found ( |Giavalisco et al 



1998 



Connolly, Sza- 



lay & Brunner 1998). Although eventually it may be possible 



pie to separate the models at the 3<r level. Because the effects 
of bias and redshift-space distortions become difficult to es- 



to avoid the uncertainties introduced by bias by measuring 
the dark matter fluctuations directly through gravitational 
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Figure 11. The power spectrum of the filtered density field p — » 
p B for B = 1.0, 1.2, 1.4 and 1.6 in the ACDM model. The y-axis 
is the ratio of A 2 (fc) for the filtered field, in real space at z = 0, 
to the BBKS prediction. All the results have been scaled by a 
constant to agree at low-fc, i.e. any linear bias has been removed. 
The thick solid line is the linear theory result, for comparison. At 
low-fc linear bias is a good approximation to these results, with 
the bias becoming more scale dependent at higher-fc. The second 
feature remains in all of these schemes. 




10 100 
K 

Figure 12. The persistence of oscillations in the 2-dimensional 
power spectrum Pqjj^K). Solid lines show the full non-linear cal- 
culation, dashed lines linear theory. The galaxy selection func- 
tions is a Gaussian centered at z = 0.3 with width 20%, which 
should be easily selectable using photometric redshifts. Note that 
the features, which become immeasurably small with an APM- 
like selection function, survive convolution with a narrow redshift 
window. The error bars approximate those expected from the 
SDSS survey. The model is the same as in Figure [id. 



lensing (Soljak 1997; Kaiser 199cj ), for the upcoming large 
surveys it will remain an issue. 

Lastly, this analysis has focused on modes with k > 
0.04 /iMpc" , since on scales larger than this several sys- 
tematic issues become important. For example, a correction 
for extinction to 10% of the true power across the survey is 
required to control systematic errors on the power spectrum 
on larger scales (for currently popular models), and this is 
quite difficult to achieve. Correlations in the photometric 
errors have similar effects. 



5.3 2D Surveys 

The discussion so far has focused exclusively on redshift sur- 
veys. However for the SDSS there will be a much larger num- 
ber of objects for which the photometry and positions will 
be known, but not the redshifts. As mentioned previously, in 
a fully 2D survey the oscillations in the power spectrum will 
become immeasurably small because of projection effects. 
But by making use of photometric redshifts, the photomet- 
ric catalogue can be used t o select redshift shells , which may 
easily be < 20% in width ( |Brunner et al. 1998| ). With such 
shells, the features in large volume 3D surveys are not signif- 
icantly washed out by projection. We shall work with shells 
of 20% width. Such relatively broad shells have the advan- 
tage of containing a large number of galaxies, thus reducing 
the shot noise. As it turns out, we will show that it is not the 
shot noise that is the limiting factor, but rather the number 
of uncorrelated measurements of the power spectrum that 
may be made. For this reason 3D surveys hold more promise 
than quasi-2D surveys for finding the baryonic features. 

A second reason to consider relatively broad shells is to 
reduce the effect of redshift space distortions. Photometric 



redshifts are a combination of true redshifts (since they are 
sensitive to the shifting of spectral breaks through photo- 
metric filters), and distance (since fluxes drop off as r~ 2 ). 
If only colors are used to compute the redshift then the er- 
rors on z axe larger, but the selection is more purely a red- 
shift measurement. Selecting galaxies on the basis of red- 
shift enhances the resul ting power s pectrum over the real 
space power spectrum (Kaiser 1987). It is straightforward 
to show, however, that if the selection function dN ga _i/dz is 
slowly varying on the scales k" 1 under consideration, one 
accurately measures the real space power spectrum. More 
exactly, for a Gaussian shell of width cr, the increase in the 
real space power spectrum (or angular correlation function) 
to leading order in ko 3> 1 is 0[f/(ka) 2 ] or smaller, where 

/(n )^ng- 6 . 

To see this, compare the fc— component of a planewave 
redshift space fluctuation <5 s (k) with the real space fluctu- 
ation <5 r (k): <5 s (k) = (1 -+- /p; 2 )5 r (k), where r ■ k = rkfi. 
The projected counts along a given line of sight are then 
SN — J ds Sn s (r) w J ds 5 s (r)p s 0(s), where 4>(s) is the 
selection function in redshift space. Using the planewave de- 
composition and retaining terms of only first order in S gives 
5N& f dkk 2 J dfj, p r <5 r (fc)(l + //i 2 ) J ds (f>(s) exp(ife^s). If 
4>{s) is a Gaussian of width cr, then J ds <j>(s) exp(ifc/is) = 
exp[— er 2 (fc/i) 2 /2]. The integration over (i then ensures 5N 
agrees with the real space projected counts to order f /(ko) 2 . 
A similar analysis starting with Limber's equation yields an 
identical result. 

On scales k ~ 0.1/iMpc -1 this correction is negligible 
for the 20% width we have chosen. If we probed larger scales, 
k — *• 0, we would eventually recover the famous Kaiser fac- 
tor (1 + 2//3 + / 2 /5), which enhances large-scale power in 
redshift space. For 20% shells this enhancement could affect 
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the first feature in P(k) near k ~ 0.02hMpc -1 , and a more 
sophisticated anaiysis wouid be needed. 

We show a 2D projection power spectrum in Fig. [lij 
The circular Sloan area has an effective radius of = 1 ra- 
dian (for a fiat sky). Its transform is 2Ji (KO)/(KO), which 
has its first zero at K = 3.83/6. We should therefore be able 
to deduce uncorrelated estimates of angular power spectra, 
provided the bins in angular wavenumber are wider than 
AK = 4. The expected fractional power errors in this case 



crp 
P 



4tt 



AKAK 



(16) 



assuming the minimal uncorrelated sampling of AK — 4. 
Note that, in 2D, the error bars decline with wavenumber 
much more slowly than in 3D, reflecting the larger num- 
ber of modes in the 3D case. This difference means that it 
is much harder to detect the baryon oscillations in the 2D 
case; for the models considered here, the characteristic sec- 
ond harmonic would never be detectable. Since shot noise 
is not a major contributor to the error bars in either 2D or 
3D, the advantage of a larger number of objects in the sam- 
ple is vastly outweighed by the reduced number of modes. 
This consideration obviously applies to other methods of 
measuring the power spectrum (e.g. lensing) where only 2D 
information is available. 



6 SUMMARY 

If the baryonic density is an appreciable fraction of the mat- 
ter density then the large-scale power spectrum is predicted 
to have features in addition to the break at the scale of 
matter-radiation equality. We have examined the observabil- 
ity of such features, specifically the oscillations on scales of 
10-100 ft" 1 Mpc which are the remnants of acoustic oscilla- 
tions in the photon-baryon fluid before last scattering. Our 
principal conclusions may be summarized as follows: 

(1) The linear theory of these processes is well un- 
derstood, and may be described as resulting in a series of 
'bumps' in power relative to a reference zero-baryon model. 

(2) The effect of largest amplitude (a factor perhaps 1.5 
in power) is a broad hump at low k. A clearer observational 
feature is the second harmonic, which is expected to lie near 
k = 0.055 Mpc -1 for most models, and which is relatively 
narrow in k space. 

(3) Non-linear effects act to wash out the oscillations at 
higher k due to mode coupling. Second order perturbation 
theory provides a qualitative description of the effect, but is 
not quantitatively accurate. 

(4) For acceptable values of the cosmological and baryon 
densities, the oscillations may be measurable with forthcom- 
ing surveys. This will best be done with three-dimensional 
surveys, rather than quasi-2D surveys based on photometric 
redshifts or lensing. The measurement of a baryon signature, 
however, will be extremely demanding. A detection of the 
broad low-fc hump, the clearest signature of the baryons, 
will require accurate measurements of the power spectrum 
for k < 0.05 ft Mpc -1 . Because of degeneracy between mod- 
els, an unambiguous detection of the second harmonic will 
require a precise determination of the amplitude of the fluc- 



tuations, as could be provided by measurements of the CMB 
anisotropy for example. 
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APPENDIX A: FIRST-ORDER 
PERTURBATIONS 

The matter and radiation power spectra used in this paper 
were calculated numerically by evolution of the coupled Ein- 
stein, fluid and Boltzmann equations in synchronous gauge 
(for details of the code see White & Scott 1996; Hu et al. 
1996). However to understand the physics behind the oscilla- 
tions it is simpler to use the Newtonian gauge. The relation 
betwee n thes e two gauges is discussed in e.g. Hu, Spergel & 
White ( [l997| ). 

We start with the fundamental equations describing the 
dynamics of a relativistic fluid. A physical interpretation 
and description of these equations can be found in Hu & 
White (1996), to which the reader is referred for more de- 
tails. The evolution of the photons and baryons in a metric 
perturbed by density fluctuations in the fcth normal mode 
is given in the Newtonian gauge as 



6 
St 



•fer-6, 

■kV b 



3$, 



for the continuity equations and 



6i = fc[9 + * 



|n 7 ] 



t(0i - Vb), 



V h = --V b + kV+T(@i-V b )/R, 



(Al) 



(A2) 



for the momentum conservation or Euler equations of the 
photons and baryons respectively. The evolution of the cold 
dark matter is the same as the baryons except for the f terms 
which describe the coupling to the photons. Overdots are 
derivatives with respect to conformal time tj = J dt/a, R = 
3pi,/4p 7 is the baryon-photon momentum density ratio, and 
f is the differential Compton optical depth. The fluctuations 
are defined as Oo = AT /T — 5 7 /4 the isotropic temperature 
perturbation, Oi the dipole moment or photon bulk velocity, 
n 7 the photon anisotropic stress perturbation, 8 b the baryon 
energy density perturbation and Vb the baryon velocity. The 
gravitational sources are the perturbation to the spatial 
curvature, and the Newtonian potential. At late times 
these are dominated by the CDM and baryons (if Qb/Qo 
is large) while at early times they are dominated by the 
relativistic species. The Einstein equations are 



(fc 2 -3iT)$ = 4irGa 2 J2 



pt5 t + 3~(pi +pi)Vi/k 



fc 2 (* + $) = SttGu 2 YtPiKi 



(A3) 
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where the sum is over particle species and the curvature 
K = -H$(l-n -n A ). 

At early times the density is high and the scattering is 
rapid compared with the travel time across a wavelength. 
Thus we may expand the momentum conservation equation 
in powers of the Compton mean free path over the wave- 
length k/f. To lowest order we obtain the tight coupling 
approximation for the evolution 



drj 3 



k 2 (l + R)^- A 3 (i + m$. (A4) 
drj 

which is a driven harmonic oscillator with natural frequency 
c~ 2 — 3(1 + R). During the tight coupling phase the am- 
plitude of the baryon perturbations cannot grow, it under- 
goes harmonic motion with an amplitude which decays as 
(1 + _R)~ 1//4 and a velocity which decays as (1 + R)~ 3 ^ 4 . 
For values of f2a consistent with standard BBN, R < 0.3 so 
this effect is not dominant. If we define the optical depth 
Ti,(rj) = fdrj'/(l + R), we find that the baryons de- 

couple from the photons when Tfc ~ 1. The oscillations 
in the baryons are frozen in at this epoch. Expanding to 
higher order in k/f one finds the oscillations are exponen- 
tially damped with characterstic scale (for more details and 
a discussion of the physics of the damping see e.g. Hu & 
White 1997) 



(A5) 



k- 2 (n) ~- fdn- ^ + 16(l + fl)/15 

where koirj) is evaluated at the peak of the visibility func- 
tion fb exp(-Tb). For Qb consistent with standard BBN this 
decoupling is after recombination. In any case, this damp- 
ing turns out not to be phenomenologically interesting, since 
non-linear effects wash out the higher peaks anyway. 

Once the photons release t heir hold on the baryons we 
can neglect the f terms in Eq. (A2). The solutions are then 
the well known growing and decaying modes of pressureless 
linear perturbation theory (e.g. in a critical density, mat- 
ter dominated universe the growing mode oc a). The den- 
sity and velocity perturbations from the tight-coupling era 
must be matched onto the growing and decaying modes for 
the pressureless components, including the CDM potentials. 
The final spectrum is the component which projects onto the 
growing mode. As is well known (e.g., Padmanabhan 1993, 
Section 8.2) at high k the growing mode is sourced primarily 
by the velocities, while at low-fc it is a mixture of density 
and velocity terms. For this reason at high k the oscilla- 
tions in A 2 (A;) are out of phase with the peaks in the CMB 
anisotropy spectrum which arise predominantly from pho- 
ton densities. In this limit peaks occur at fer* = (2j + l)n/2 
where j = 0, 1, 2, • • • and r* is the sound horizon at decou- 
pling: r* = J c s drj. In detail this differs from the sound 
horizon at recombination (which controls the position of the 
CMB peaks) but as can be seen in Fig. [j] the two horizons 
are comparable. 

In Fig. [j] we see that the baryons induce both oscilla- 
tions and amplitude suppression in T(k). Let us consider the 
latter first. Modes which enter the horizon before matter- 
radiation equality can grow only logarithmically at best. 
For the baryons, modes which enter the horizon after equal- 
ity but before decoupling oscillate with decaying amplitude. 
They therefore do not contribute fully to the gravitational 
potentials in Eq. (A3). For the CDM the modes which enter 



the horizon after equality but before decoupling grow as if in 
a universe with density » fio — rather than ilo- Hence for 
the period between equality and decoupling all modes which 
are inside the horizon have their growth suppressed relative 
to the Ab *C Oo case. This damping of power changes the 
shape of P(k) near its peak, with the largest effect occuring 
in models where equality and decoupling are most separated 
(high fish 2 and high Q,oh 2 ). 

The acoustic oscillat ions in the baryons at decoupling 
(Hu & Sugiyama 1996) are then superposed upon the 
smooth power spectrum. The amplitude of the oscillations 
depends on both the driving force ($ and ^) and on R. 
Since the potentials decay in the radiation dominated epoch, 
larger oscillations come from higher fish 2 and lower Q,oh 2 
(see e.g. discussion of the potential envelope in Hu & White 
1997). Specific examples are shown in Fig. ]4. 



APPENDIX B: PM CODE 

To investigate the effects of non-linearity on the persistence 
of the features in the power spectrum we used a PM code 
to calculate the non-linear power spectrum. Since it has not 
been discussed before, we give some details of our implimen- 
tation of the code here. All of our results were also checked 
by running P 3 M simulations as discussed in the text. Above 
the grid scale the agreement was excellent. 

We used 128 particles with the gravitational forces 
computed on a 128 3 grid. All of the runs were started at 
1 + z — 20, and evolved using equations of motion in which 
the 'time' coordinate is In a, where a is the scale factor. The 
initial particle positions were displaced from a random posi- 
tion within a cell (with cells uniformly filling the box) using 
the Zel'dovich approximation. Each realization of the power 
spectrum was chosen to have random phases and amplitude 
drawn from an exponential distribution. 

To calculate the forces on the particles we assigned them 
to a grid using the cloud-in-cell algorithm. On large scales 
the power spectrum recovered is independent of the charge 
assignment scheme, which we explicitly checked. The forces 
were computed using FFT techniques with a force kernel of 
k/fc 2 to compute V$. While it is computationally faster to 
compute $ directly using a l/k 2 kernel and calculate V<3? 
by differencing, we found that this gave worse performance 
on small scales than calculating V<3? for each direction using 
k/fc 2 . The better accuracy at small scales in this method 
more than makes up for the extra computing time (com- 
pared to running higher resolution). The time step was dy- 
namically chosen as a small fraction of the inverse square 
root of the maximum acceleration, with an upper limit of 
Aa/a — 4% per step, where a is the cosmological expansion 
factor. This resulted in a final particle position error of less 
than 0.1% of the box size. 

At selected time steps the power spectrum was calcu- 
lated. The density field 5k was computed by assigning all 
particles to the nearest node in a 128 3 grid and perform- 
ing an FFT. The resulting power, |<5fc| 2 , was corrected for 
shot-noise and binning onto the grid to obtain an estimate 
for P{k). The redshift space power spectrum was computed 
by adding one component of the velocity (in units of the 
Hubble constant times the length of the box) to the particle 
positions before assigning them to the grid. 
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APPENDIX C: SCALE-DEPENDENT BIAS 



Let us imagine that we have two classes of objects, whose 
probability of being included in the sample at a position 
R, is </> a (R.) with a = 1, 2. (Luminosity is only one possible 
criterion for inclusion in a catalog; surface brightness could 
be another.) A straightforward calculation gives 



4>li4>lj^ll + (01i02j + 4>2i<t>ij)£,12 + 4>2i4>2j&2 



(CI) 



{4>u + 4>2i){(f>ij + 4>2j) 

where i, j label positions in the survey and we have used the 
shorthand (fiij = 0i(Rj). If we further suppose that ^n(r-) = 
bU(r), £ 22 (r) = &|£(r) and £ 12 (r) = &i& 2 £(r), where £(r) 
represents the correlation function for points at a separation 
r drawn from the underlying field, of which both populations 
are biased tracers. We find 



b 2 (R, 



= N- 



(6i0i + &202)i(6l01 + fr2<fe)j 
'01 + <f>2)i{4>l + <t>2)j 



(C2) 



(C3) 



where i? is the depth of the survey, the averaging is over the 
survey area, the sum is over all positions i,j in the survey 
which are separated by a distance r and iVp r is the number 
of such pairs of positions. The resulting power spectrum will 
be 



A 2 (R,k) 



(2tt) 3 



d 3 g& 2 (fl,|k-q|)A§(ii, g ) 



(C4) 



where b 2 (R, k) is the Fourier transform of b 2 (R, r) and A§(fc) 
is the power spectrum assuming no bias. 

If the selection functions vary smoothly throughout the 
survey volume, the inferred bias parameter will depend only 
on the survey depth R as the ratio <j>i(R)/<j)2(R) varies with 
R. This results in an overall shift with depth of the ampli- 
tude of the measured power spectrum, which would need 
to be distinguished from evolution in the clustering. It is 
also possible, however, to introduce a feature into the power 
spectrum at a given depth if the selection function of one of 
the populations changes rapidly in the survey volume. This 
could occur in a redshift slice if the selection function of one 
population either was rapidly declining in the slice or oc- 
cupied a very narrow range in redshift space. For galaxies 
at a fixed mean redshift in a very narrow redshift shell, the 
population composition that dominates the contribution to 
the power spectrum at a given scale will differ for different 
scales. Depending on the shapes and peaks of the selection 
functions, and the position of the shell, this can result in 
either enhancements or reductions in the power spectrum 
with increasing scale. 



We illustrate the effect in Figure CI. The survey vol- 
ume is taken to be cone with solid angle -k sr, as for the 
SDSS. The selection functions for both populations are as- 
sumed to be Gaussians with mean redshift Zi and width oi. 
One population is taken to comprise 75% of the galaxies 
(integrated over all redshifts), with z\ = 0.2, o\ — 0.1, and 
&i = 2. The selection function for the remaining population 
is very narrowly peaked at z% = 0.3 with 02 = 0.003. It 
is taken to be unbiased (62 = 1). The underlying matter 
power spectrum is assumed to be BBKS CDM. For a shell 
centered at z = 0.3, an oscillating scale-dependent bias is 
produced, which could easily be mistaken for an acoustic os- 




0.1 

k (h Mpc" 1 ) 



Figure CI. Scale— dependent bias is introduced by mixing a 
dominant population having a broad selection function peaking 
at Z\ = 0.2 and bias bi = 2, with a minority unbiased (62 = 1) 
population having a selection function that is sharply peaked 
at Z2 = 0.3. The resulting distortion of the assumed underly- 
ing BBKS power spectrum is shown for redshift shells positioned 
at z = 0.3 (solid), z = 0.295 (long-dashed), z = 0.293 (short- 
dashed), and z = 0.292 (dotted). The features could easily be 
mistaken for acoustic oscillations at z = 0.3, and would suppress 
the oscillations at z = 0.292. They would, however, be averaged 
away in a shell of full width Az = 0.1. 



dilation. The effect, however, diminishes sharply as the shell 
is moved away from z — z 2 , until by 2 = 0.25 the distinc- 
tions are reduced to less that 1%. This means that for a shell 
of full width Az — 0.1, the oscillations would be completely 
masked. 

The strength of the effect is entirely a consequence of 
the narrowness of the selection function for the second pop- 
ulation, which we consider to be unphysical. Moreover, even 
if the effect were detected, its extreme sensitivity to redshift 
would distinguish it from the acoustic oscillations, which are 
not expected to change with redshift except by the much 
more gradual growth due to gravitational instability. Typi- 
cally, the effect on the power spectrum is much smaller than 
that shown here, even in the presence of a substantial scale- 
dependent bias. For instance, if the width of the second pop- 
ulation above were increased to 02 = 0.05, with everything 
else held unchanged, at z — 0.3 the distortions of the power 
spectrum would be only 0.1%, even though the mean bias 
factor changes by 3% between 20 and 200 ft -1 Mpc. 



APPENDIX D: GALAXY NUMBER COUNTS 

The number of galaxies per unit redshift dN/dz to a limiting 
magnitude mn m in a survey of area Q s sr is related to the 
proper luminosity function cj>{M, z) of the galaxies by 
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dN , . _ l2 , , dL 

-5— (< mum) = il s d A (z) — 
dz dz 



M max 



dM 4>(M,z) 



where 

M max : 



mii„ 



51og 10 d h (z) - 25 -K{z)- A, 



(Dl) 



(D2) 



di,{z) is the luminosity distance (in Mpc) to a source at 
redshift z, c£a(z) = dh(z) / (1 + z) 2 is the angular distance to 
the source, and d£ p is the proper differential line element. 
The terms K(z) and A account for the K-correction and 



extinction, respectively, lor A = U, 



*,(*) = {zVo - (2 - Q ) [(1 + fi z) 1/2 - l] } , (D3) 



and 



^ = ^(l + ,)-(l + n o«)-^. 



(D4) 



When A / closed form expressions do not exist and the 
defining integrals need to be evaluated numerically. 

For a Schechter luminosity function, log 10 (L/L*) = 
0.4(M* - M) and 

<f>(M,z) = (0.4 log 10) {1 + zf <j>*(z) 
x(L/L,) 1+a exp[-L/L»], 



(D5) 



where the density expansion rate (1 + z) 3 has been factored 
out so that any residual redshift dependence in <f>* (z) is due 
to galactic evolution. The integral over M is now analytic. 
Defining x min by 2.5 log 10 x min = M* - M max , 



dM <j>(M,z) = (l + z) 3 <f>*(z)r(l + a,x min ), (D6) 



where V(a, x) is an incomplete gamma function ( Gradshteyn 
& Ryzhik 1980). For — 2 < a < — 1, it may be evaluated 
using the recursion relation, 



r(i + a ,x) = — )— [r(2 + , 

1 + a 1 



(D7) 



We have used these estimates to obtain the number of 



galaxies in redshift shells in Section 5.S and hence the shot- 
noise contribution to the error bar in Fig. [l^ . 
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